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Abstract 

We present the results of a partial upgrade to the Monte Carlo event generator TAUOLA using Resonance Chiral 
Theory for the two and three meson final states. These modes account for 88% of total hadronic width of the tau 
meson. The first results of the model parameters have been obtained using Preliminary BaBar data for Tin mode. 
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1. Introduction 

The tau lepton is the heaviest of the leptons, and pro- 
vides a unique opportunity to study low energy QCD 
and the mechanism of hadronization. With the shut- 
down of both the B-factory collaborations, BaBar |[T] 
and Belle [2J, this is a critical time to make the results 
from more than a decade of experimental research avail- 
able in a useful manner to high energy physics commu- 
nity before the opportunity is lost. To accomplish this, 
interaction between experimentalist and theorists is re- 
quired to determine the optimal way for comparing ex- 
perimental data with the theoretical prediction. The ap- 
proach used by this work is to compare theory to data 
by using a Monte Carlo event generator to simulate the 
theoretical models. 

The Monte Carlo event generator TAUOLA is a long 
term project that started in the 90's with the publication 
The generator simulates more than twenty modes, 
including both the leptonic and hadronic modes. Mod- 
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eling the hadronic decay modes involves matrix ele- 
ments that convey the hadronization of the vector and 
axial-vector currents. At present there is no determi- 
nation from first principles for those matrix elements 
since they involve strong interaction effects in the non- 
perturbative regime. Therefore, one has to rely on mod- 
els that parameterize the form-factors originating from 
the hadronization. The hadronic currents implemented 
in the generator TAUOLA Q were based on theoretical 
results presented in |4|. 

In these models, the hadronic form-factors are writ- 
ten as a weighted sum of products of Breit-Wigner func- 
tions. This approximation, as it is demonstrated in [5|, is 
not able to reproduce the next-to-leading-order ;^'PT re- 
sults. Later the experimental collaborations, both Cleo 
||6] Ul and Aleph [8 |, introduced improvements based 
on the result of their data analysis and in some cases 
spoiling the theoretical constraints [^J As an alternative, 
it was proposed to apply the methods of the Resonance 



In fact, hadronic currents of Cleo and Aleph versions spoil some 
theoretical prediction. For example, to reproduce t" — > (KKnYvj 
data the Cleo collaboration Q introduced two ad-hoc parameters and. 
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Chiral Theory Q. This approach is a better theoreti- 
cally founded than the approach of weighted products 
of Breit-Wigner functions. However, it requires com- 
parison both BaBar and Belle experimental data with 
the results from theoretical models which have to be 
implemented into TAUOLA. In the next sections, we 
will present the result on the implementation the RChT 
hadronic currents for n^jfi, K K^\ {KnY , {nrniY and 
{KKn)- modes in TAUOLA, for details see IfTOlfTTl . 

2. Hadronic currents for two and three meson decay 
modes 

For any two meson r decay channel, the most gen- 
eral form of hadronic currents (y^), which is compatible 
with Lorentz invariance, is written as 

J'' = N\{p,-p2-—{pi+P2)fF''{s) (1) 

L S 

+ —(i.Pi+P2YF\s)\ 
s J 

where pi and p2 are the momenta of hadrons, A12 - 
m\- rn^, s - {p\ + pif"- For a final state of three pseu- 
doscalarsis 

Ji" = N[T>i[c,{p2 - p.YFi + C2(p3 - (2) 

where pi, p2 and ^3 are the momenta, T^^y = - 
Itilv/'f' denotes the transverse projector, - (pi + 
Pi + PiY is the momentum of the hadronic system and 
F is the pion decay constant in the chiral limit. Among 
the form factors F\, F2, f 3, corresponding to the axial- 
vector part of the hadronic currents, only two are model- 
independent. For our convenience we keep all of them 
in Eq. (|2|. 

The model-dependence in Eqs. ([T]i and (j2]l is in- 
cluded in the hadronic form-factors {Fy, Fs expressed 
in terms of F,, i = 1...5). For three meson decay modes, 
the hadronic form-factors calculated within Rj^fT can be 
written as 

Fj^F] + Ff + f (3) 

where F^j is the chiral contribution, F^ is the one res- 
onance contribution and is the double-resonance 
part. The explicit form of the functions Fi for 3n and 
KKn modes can be found in ifTOl . Section 2, and in 
ll5l fT2]| . The corresponding theoretical form-factors are 



as a result, the Wess-Zumino part does not reproduce the QCD nor- 
maUzation. 



obtained in the isospin limit {m,^ — {2m„+ + m„o)/3, 
niK - {niK* + nij^o)/!), except for the two pion and two 
kaons modes. 

3. Numerical results. Fit of the three pion mode to 
BaBar data 

To avoid problems with multi-dimensional integra- 
tion of the fli -meson propagator, which is rapidly- 
changing as a function of its arguments, we first tabu- 
lated the Faiiq^) of Eq. (38) of IfTOll . Then we use linear 
interpolation to obtain the value of the a; width at the 
required q^. 

To check the numerical stability of the generator and 
the multiple numerical integration, the following tests 
have been done: 

• For every channel one dimensional spectrum 
dr/dq^ (for the three meson decay modes) and 
dYjds (for the two meson decay modes) pro- 
duced by the generator has been compared with 
the semi-analytical results. For the three meson de- 
cay modes the Gauss integration method has been 
applied to integrate the analytical results for the 
hadronic currents. The results agree within statis- 
tical errors except for the first and last bins. 

• The spectrum obtained by the Gauss integration 
has been compared with the linear interpolated 
spectrum from the neighboring points and demon- 
strated that the fluctuations due to numerical prob- 
lems of integration are absent; 

• The total rate for every channel obtained from a 
Monte Carlo run has been compared with the semi- 
analytical method and required to agree within sta- 
tistical error. The results of comparison are pre- 
sented in TablT] 

The tests have been performed in the isospin limit of 
meson masses, except for the two pion and two kaons 
modes. For more details about these tests see FTOl. 
Triple Gaussian integration is used for the analytical cal- 
culation and double Gaussian integration is used for the 
current calculation that enters the matrix elements of the 
Monte Carlo generation. Thus, a pre-tabulation of the 
cii width, Ya^iq^), is a convenient tool to significantly 
increase the generation speed. 

To obtain the proper kinematic configurations, the 
differences between neutral and charged pion and kaon 
masses were taken into account, more specifically, the 
physical values were chosen in the phase space genera- 
tion. On the other hand, this choice breaks constraints 
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Channel 


Width, [GeV] 




PDG 






Equal masses 






Phase space 
with masses 






n-n' 


(5.778 ± 0.35%) 


10 


-li 


(5.2283 ± 0.005%) ■ 


10- 


13 


(5.2441 ± 0.005%) 


10- 


-13 


K-n° 


(9.72 ± 3.5% ) ■ 


10- 


15 


(8.3981 ±0.005%) ■ 


10- 


15 


(8.5810 ±0.005%) 


10- 


-15 


n-K^ 


(1.9 +5% )■ 


10- 


14 


(1.6798 ± 0.006%) ■ 


10- 


14 


(1.6512 ±0.006%) 


10- 


-14 




(3.60 + 10% ) ■ 


10- 


15 


(2.6502 ± 0.007%) ■ 


10- 


15 


(2.6502 ± 0.008%) 


10- 


-15 




(2.11 ±0.8% ) 


10- 


-13 


(2.1013 ±0.016%) ■ 


10- 


13 


(2.0800 ± 0.017%) 


10- 


-13 




(2.10 ± 1.2% ) 


10- 


-13 


(2.1013 ±0.016%) ■ 


10- 


13 


(2.1256 ±0.017%) 


10- 


-13 


K-7T-K+ 


(3.17 ±4% )■ 


10- 


15 


(3.7379 ± 0.024%) ■ 


10- 


15 


(3.8460 ± 0.024%) 


10- 


-15 




(3.9 +24% )■ 


10- 


15 


(3.7385 ± 0.024%) ■ 


10- 


15 


(3.5917 ±0.024%) 


10- 


-15 




(3.60 ± 12.6% ) 


■ 10 


-15 


(2.7367 ± 0.025%) ■ 


10- 


15 


(2.7711 ±0.024%) 


10- 


-15 



Table 1: The r decay partial widths. The PDG value 1 13], the 2nd column, is compared with numerical results from TAUOLA with the R^T cuiTcnts 
in the 3rd column (obtained in the isospin limit for meson masses) and the 4th column (using physical masses for the pahse space). 



resulting from isospin symmetry in potentially uncon- 
trolled way. This is why we collect the numerical re- 
sults from the Monte Carlo calculation, shown in Table 
[T] where the partial widths from Particle Data Group 
compilation 1 1 3 1 are compared with our results obtained 
with isospin-averaged pseudoscalar masses and with the 
physical ones. The model parameters, more specifically 
the masses of the resonances and the coupling constants, 
were fitted to Aleph data (14], requiring correct high- 
energy behavior of the related form factors. For details 
see Appendix C in IfTOl . From Table [T] one can see that 
the difference between the TAUOLA results for the r 
decay partial width and the PDG ones is 1.5% - 17% 
depending on the mode. As expected, the agreement is 
not good because only minimal attempts to adjusting the 
model parameters have been applied for the comparison 
with BaBar and Belle data. 



Currently, only the differential spectrum of the two 
pion modes f2] and three pseudoscalar modes (15] are 
published. We begin with a fit to the spectrum of the 
n^n^n^ mode. The result is presented in Table |2] and 
Fig[3] The fit was done taking into account only the 
dominant S-wave mechanism (we follow the determi- 
nation done in lfT6l ). As suggested in lfT6l the discrep- 
ancy in the low mass region could be described adding a 
contribution from a scalar particle, f -wave mechanism. 
We expect that inclusion of the lowest-lying scalar res- 
onance (T7] will improve the value ofx^- The values in 
the 5th row of Table |2] are only the preliminary results. 
They do not necessarily correspond to the minimum of 

Indf of the final fit. Work is in progress. Some tech- 
nical aspects of the fitting strategy is given in talk by Z. 
WaslfTSl. 



4. Tau physcis at LHC 

From the perspective of high energy experiments, 
such as those at LHC, a good understanding of tau lep- 
tons properties contributes important ingredients of new 
physics signatures. With the discovery of a new parti- 
cle around the mass of 125-126 GeV |19|, tau decays 
are an important decay mode for determining if this 
is the Standard Model Higgs. This is especially per- 
tinent since CMS has reported a deficit in the number 
of fermion decays from the new particle relative to the 
Standard Model Higgs Prediction. 

At LHC, at the moment, tau decays are only used for 
identification and are not used to study their dynamic. 
However, the dynamics of tau decays are important for 
both modeling the decays and therefore the reconstruc- 
tion, identification and for measuring the polarization 
of tau decays. As a result, an upgrade to the TAUOLA 
based on the BaBar and Belle measurements of tau de- 
cays is of some interest for systematic errors evaluation 
for LHC measurements. 

5. Conclusion 

The theoretical results for the hadronic currents of 
two and three pseudoscalar modes, namely, nn, Kn, 
KK, nnn and KKn, in the framework of R;^'T have 
been implemented in TAUOLA. These modes, together 
with the one-meson decay modes, represent more than 
88% of the hadronic width of the tau lepton. Rj^^T is 
a more controlled QCD-based model than the usually 
used Breit-Wigner parameterization. However, before 
making conclusion about validity of the model the the- 
oretical results have to be confronted with the experi- 
mental data. This can be achieved by fitting the model 
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Data - BaBar Preliminary Tau 2012 
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Figure 1: Invariant n^n n (left) and n^n (center) and ir n (riglit) mass distributions. Lighter grey liistograms are for R^T parametrization, 
darker grey is for CLEO one, points are Babar data. 





Mp, 




M„ 


F 


Fv 


Fa 




Min. 


1.44 


0.32 


1.00 


0.0920 


0.12 


0.1 


-0.36 


Max 


1.48 


0.39 


1.24 


0.0924 


0.24 


0.2 


-0.18 


Default 


1.453 


0.4 


1.12 


0.0924 


0.18 


0.149 


-0.25 


Fit, x^/ndf = 2262.12/132 


1.4302 


0.376061 


1.21706 


0.092318 


0.121938 


0.11291 


-0.208811 



Table 2: Numerical values of the R^T parameters fitted to BaBar data for three pion mode 1151 



to data. As a consequence of comparisons with the data, 
some of the theoretical assumptions may need to be re- 
considered too. Now that the technical work on current 
installation is complete, the work on fitting the data is in 
progress in collaboration with theoreticians and experi- 
mentalists. Therefore, we consider this work as a step 
towards a theoretically rigorous description of hadronic 
tau decay data. 

As soon as the theoretical results are able to repro- 
duce Babar and Belle data for the most important for 
LHC decays into two and three pions the TAUOLA- 
RChT version will be installed into LHC environment 
as it is described in ["201. 
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